A lattice study of the masses of singlet ++ mesons 
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We compute the masses of the flavour singlet ++ mesons using (n/ = 2) unquenched lattice QCD 
with the fwasaki and Wilson gauge actions. Both fermionic and glueball interpolating operators are 
used to create the states. The mass of the lightest ++ meson is suppressed relative to the mass 
of the ++ glueball in quenched QCD at an equivalent lattice spacing. We discuss two possible 
physical reasons for this. 
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I. INTRODUCTION 



The interpretation of the experimental fo mesons in 
terms of fundamental quark and glue fields is still not set- 
tled QLulbiQlLAuI- Quenched lattice QCD predicts that 
the mass of the scalar {J PC = ++ ) glueball is around 
1.6 GeV Hi, SGI. Hence, attention has focused on 
finding evidence for a glueball component (mixing with 
quark-antiquark) in the physical / (1370), /o(1500), and 
/o(1710) mesons. The quark model predicts that there 
are two fo mesons in this mass regime, so the existence of 
three mesons is suggestive of the presence of additional 
degrees of freedom such as the elusive ++ glueball. 

The experimental fo spectrum contains more puzzles. 
The experimental data for the /o(1370) still seems con- 
troversial \M . There is also a new state /o(1790) reported 
by BES [jjj . It is not clear how this effects the "stan- 
dard" mixing scenario for the ++ mesons, that consider 
only three mesons. The interpretation of the /o(980) me- 
son which is close to the KK threshold, in terms of quark 
and antiquarks, is also uncertain. Even extracting the 
mass and width of the fo (400-1200) state is still con- 
troversial 01 0, but progress seems to have been made 
recently [l^, uM ■ ^ may be difficult for lattice calcula- 
tions to explore the /o(400-1200), because it has such a 
large width. Although there are attempts to study this 
state |3 E3 on the lattice. 

The glueball spectrum in pure gauge theory is the- 
oretically well defined, because the glueball operators 
do not mix with fermionic tpip operators. The work of 
many authors has shown that the lightest ++ state is at 
1640(40)MeV (plus errors of around 10% in setting the 
scale) 

in pure SU(3) Yang-Mills theory. In un- 
quenched calculations there are new complications. The 
glueball and fermionic ++ states can mix as illustrated 
in figure ^ 111 fact, it no longer makes sense to talk 
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FIG. 1: Mixing as unquenching occurs. 



about glueball states in dynamical QCD calculations (al- 
though there are glueball interpolating operators). There 
are only flavour singlet ++ mesons that we sometimes 
denote by FS. Also the ++ states can decay into meson 
pairs, so the decay width of hadron may play an impor- 
tant role in the dynamics. 

Unlike in quenched QCD which has a degenerate 
isoscalar and isovector scalar meson (made from qq), the 
isoscalar state mixes with the glueball, being reduced in 
mass if it is lighter than the glueball. In detail, this mix- 
ing will be a function of the lattice parameters and sea 
quark masses. These will have to be extrapolated to their 
physical values to get the physical mixing. 

There are predictions [TtIIT^ I that the width for glue- 
ball decay to two mesons is large relative to its mass. 
Although one exploratory lattice calculation [l^ found 
that the decay width of the ++ was 108(29) MeV in the 
quenched approximation. 

Indeed, the standard formalism for determining the 
masses of hadrons on a lattice may be inappropriate in 
the presence of open decay channels and special tech- 
niques may be needed ^ El E2- Tne MILC collab- 
oration reported problems in extracting the masses of 
non-singlet ++ |2^ and 1 h [24| mesons, that they at- 
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tributed to open decay channels. The techniques to study 
decay widths using lattice QCD have recently been re- 
viewed [iil ■ All the above discussions indicate that it is 
essential to study the singlet ++ mass spectrum with 
dynamical fermions. 

The basic formalism for the mixing of pure glue and qq 
states was described in previous work [26| from UKQCD. 
That paper obtained surprisingly light values for the 
flavour-singlet scalar meson, but at relatively coarse lat- 
tice spacing. Here we explore this issue using a smaller 
lattice spacing. We also determine the spectrum using 
lattices with a different gauge action since the lattice ar- 
tifacts (for instance big order a 2 corrections to the scalar 
mass) will then be different. 

In detail, we use dynamical gauge configurations gen- 
erated by UKQCD using fully non-perturbatively im- 
proved clover fermions [27|. The lattice spacing is finer, 
a ~ 0.1 fm, than in our previous study [2(J. Another 
technical improvement is that we use fermionic scalar 
mesons at nonzero momentum. This helps to get a bet- 
ter signal for the larger volumes. This was regularly used 
for pure glue operators 28] , but is not routinely used for 
fermionic operators. The glueball spectrum using only 
glue operators has been presented in [53, H^]. The in- 
clusion of qq operators in this study allows us to address 
mixing issues. UKQCD has recently published a lattice 
study of the non-singlet ++ mesons |29j . 

We also use gauge configurations from the CP-PACS 
collaboration [3(|. These configurations were generated 
with an improved gauge action. This can help to address 
some of the potential issues with using the clover action in 
combination with the Wilson pure gauge action [M], [3^] . 
In particular the lattice spacing dependence of the mass 
of the ++ states is known to be large with the Wil- 
son pure gauge action ||| 13, |U ■ There have also 
been claims that the unqucnchcd calculations that use 
the clover fermion action, with the Wilson single plaque- 
tte action, are affected by bu lk p hase transition in the 
adjoint plane 0, 113 H3, H3, |2a 1^3 . One conjectured 
consequence of the phase transition was the suppression 
of the ++ glueball masses obtained by Hart and Te- 
per [H from Nj = 2 unquenched QCD. The JLQCD 
collaboration |3y| found the effect of the adjoint phase 
transition was reduced by the use of improved gauge ac- 
tions, such as the Iwasaki action. Hence the CP-PACS 
data will be an important cross check on results. 



II. PARAMETERS OF THE LATTICE 
CALCULATION 

We use gau ge c onfigurations from the UKQCD j2?J 
and CP-PACS 30] collaborations. For the gauge config- 
urations from UKQCD, the non-perturbatively improved 
clover action was used to generate the unquenched gauge 
configurations, with a clover coefficient csw = 2.0171. 
The (3 value was 5.2 and the lattice volume was 16 3 x 
32. The sea quark k values of 0.1350 and 0.1355 were 



used for the singlet correlators. The hadron spectrum, 
potential, and some glueball estimates from this data set 
have been presented in [23, H3] • 

We use two sets of data with Nt = 2 from the CP- 
PACS collaboration at (3 = 1.95 30]. The tadpole im- 
proved clover action with clover coefficient of csw — 1-53. 
The Iwasaki renormalised group improved gauge action 
was used. The lattice size is 16 3 x 32 and we use va- 
lence quark masses equal to those of the sea quarks. The 
lattice details are summarised in Table Q] 



Code 


■^V gauge 




ro/a 


amps 


C390 


648 


0.1390 


2.651 


0.729 


C410 


490 


0.1410 


3.014 


0.427 


U350 


144 


0.1350 


4.75 


0.405 


U355 


416 


0.1355 


5.04 


0.294 



TABLE I: Simulation details for CP-PACS and UKQCD data 
sets. 



The methods we use to compute the disconnected dia- 
grams have been described in earlier publications [2^.l3^ |. 
Essentially we use random Z2 volume sources to estimate 
the bubble diagrams, with the variance reduction tech- 
nique described in |2(|. We have also published data for 
the singlet pseudoscalar channel |38] from some of this 
data set. We have presented preliminary results for a 
subset of the data used in this paper in 39] . 

III. FIT METHODS 

To extract a good signal from lattice QCD calculations, 
it has been found to be essential to use a variational basis 
of correlators. 

a ii (t) = <o|o,(t)to i (o)|o) (i) 

In terms of the path integral the matrix dj is 

Cij(t) = ^ JdU J dip I ' dQe-S'-SoOiitfOjiO) (2) 

where Sf is the clover fermion action (lattice approxi- 
mation to the Dirac Lagrangian) and Sg is Wilson gauge 
action (lattice approximation to the gauge action). The 
full details of the Lagrangians are described in our earlier 
work on the light hadron spectrum [53 . In the fermion 
sector we use fuzzed and local operators as basis states 
(Oi(t)). The fuzzing method is described in the spectrum 
study [23 of UKQCD. In the pure glue sector, two types 
of smeared glueball operators j2^ are included. 

We use factorising (or variational) fits to extract the 
masses and amplitudes. 

M 

C i3 (t)=Y^Cc^e- E - t (3) 

m=0 
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TABLE II: Fit results for flavour singlet scalar mesons. These 
are fits to four by four matrices of correlations to the f-range 
shown with momentum (p — 2nn/L. Fits to momentum zero 
are to the vacuum subtracted correlators here. For these fits 
involving disconnected contributions, we expect a x 2 /dof close 
to 1.0 since the absolute error stays constant with increasing 
t. 



Code n 


Region x 2 /dof 


E 


£i 


E 2 


U350 000 


2 - 


6 


24/40 


0.64(3) 


1.17(25) 




U350 100 


2 - 


6 


24/35 


0.69(4) 


0.86(7) 


1.23(29) 


U355 000 


2 - 


6 


24/40 


0.510(14) 


1.06(12) 




U355 100 


2 - 


6 


34/35 


0.62(3) 


1.01(5) 


1.47(10) 


U355 110 


2 - 


6 


36/35 


0.75(2) 


1.05(3) 


1.63(9) 


C390 000 


2 - 


6 


45/45 


0.98(4) 






C410 000 


2 - 


6 


19/40 


0.67(2) 


1.56(40) 





The scalar channel (0 ++ ) at zero momentum has a vac- 
uum contribution which we account for either by fixing 
Eq = [40j or by evaluating vacuum subtracted corre- 
lations. M is the number of states in the fit which can 
be 1 to 3. For our multi-state fits, we quote the masses 
for all states. However we regard the highest state as a 
"truncation error" , and this mass may not correspond to 
a physical state. 

Any state with the same quantum numbers as the op- 
erator Oj will couple to that channel. If the amplitude is 
small then it will be difficult to extract the mass of the 
state. We feel it is better to use as many different inter- 
polating operators in constructing the correlator matrix. 
If the singlet ++ states couple to both pure glue and 
fermionic operators, then it must be better to use addi- 
tional basis states, provided they are not too noisy and 
are sufficiently independent of the other states. This can 
help to stablize the multi-state fits. Hence our best re- 
sults are from fits from order 4 smearing matrices that 
include both glueball and fermionic operators, each with 
two spatial sizes. 

The Iwasaki action has ghost states that contribute 
to the correlators at small time distances |33l l4Tl | . This 
issue has recently been studied by Necco [23 m quenched 
QCD for a variety of improved gauge actions. Here we 
shall restrict our fits to t > a which should reduce these 
potential problems. 



IV. RESULTS 

In figure|2]we show an effective mass plot for the U355 
data at zero momentum together with a fit. This illus- 
trates the quality of the data and also the advantage of 
having many different correlations to fit simultaneously. 

In Table |Hj we collect our fit results for the flavour 
singlet scalar. Results for the flavour non-singlet were 
presented previously [2^. These results use the 4 by 4 
matrix of correlators that includes pure glue and pure 
fermionic operators in the same fit. For connected cor- 
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FIG. 2: Singlet scalar meson effective mass plot from 
fermionic and gluonic operators from U355 with momentum 
zero. 

relators we use the conventional method with sources at 
the origin (with 4 time sources for U355). For discon- 
nected fermionic loops we use the method of j2|| with 
100 stochastic samples. For "glueball" operators we used 
square fuzzed Wilson loops of differing sizes, selected to 
have good overlap with the scalar glueball 28] . We mea- 
sure glueball operators every 10 gauge configurations and 
fermionic disconnected operators every 10 (40 for U350 
and 20 for U355). The data are binned into groups of 40 
trajectories to avoid the impact of any autocorrelation 
on the error analysis. 

We use correlated fits |43| to ascertain the minimum t- 
value that gives an acceptable fit. We then quote results 
from an uncorrelated fit to that range. We find similar 
results for the two lightest singlet states at non-zero mo- 
mentum from 3 state fits to t ranges of 2-6 to 2-10 and 
from 2 state fits for t ranges of 3-6 to 3-10. Thus the 
systematic error on the ground state energy from vary- 
ing the fits seems to be smaller than the statistical error 
which is quoted in the table. 

For the FS channel, the non-zero momentum results 
are more stable since they do not involve a vacuum sub- 
traction and also they involve additional statistics from 
constructing the correlator from momenta in each spatial 
direction. We expect that 

E 2 = m 2 + p 2 (4) 

where p = 2m: /L with integer vector n where we explore 
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FIG. 3: Scalar meson masses (for U355) versus momentum for 
both flavour singlet (/o) and non-singlet (oq). The straight 
lines show the expected slope for a relativistic mass-energy 
relationship. The lowest mr thresholds are also illustrated. 



n.n = 0, 1, 2 here. As shown in fig. 01 we find consistency 
(U355 is plotted since for that case we have 3 momentum 
values) within our statistical errors. The mass value is 
then the intercept and is consistent with the results in 
table EH 

Using glueball operators at zero and nonzero mo- 
mentum, UKQCD |2|| previously found aM 0++ = 
0.628(30) and aM 0++ = 0.626(41) at k = 0.1355 and 
at k — 0.1350 respectively. The masses were estimated 
using effective masses after finding the optimal basis us- 
ing a variational technique |28| . This previous analysis 
claimed that the ++ state was lighter in Nf = 2 with 
the mass ratio from dynamical to quenched lattices of 
0.85(3). 

Using equation 0] with the data in table [HI we ob- 
tain am 1 (0++) = 0.51(2) for U355 and 0.60(4) for U355. 
Hence, doing a combined fit to glueball and qq ++ op- 
erators has produced additional suppression of the ++ 
mass with the inclusion of dynamical fermions. It is only 
for U355 that we have a reliable estimate of the first ex- 
cited state, obtaining am2(0 ++ ) = 0.92(6). For U355, 
we find that the fit coefficients of our two lightest states 
have the structure expected from a maximal mixing of a 
gluonic state and a fermionic state. The previous analy- 
sis m using only gluonic operators appears to have 
been getting close to the weighted mean of the first and 



FIG. 4: Compilation of FS ++ masses versus a 2 . 



second masses which is consistent with the maximal mix- 
ing we find here. 

In figure 0] the ++ masses (from the four basis fit) are 
plotted with the continuum limit of the quenched data 
(with standard Wilson glue). The masses are in units of 
7"o, that is convenient unit for lattice studies. To convert 
to MeV the reader can use 1 /r — 400 MeV. 

Our study using the CP-PACS lattices with Iwasaki 
glue has revealed somewhat similar results to those ob- 
tained previously using Wilson glue [2(j at similar lattice 
spacings. Our results (U355 and U350) with finer lattice 
spacings also illustrate the suppression of the mass of the 
FS ++ state relative to the quenched ++ glueball mass 
at a similar lattice spacing. Thus the feature found previ- 
ously 26] of surprisingly light /o mass is consistent with 
the new results at a finer lattice spacing and with an im- 
proved gauge action at the coarse action. We discuss this 
further in section IVTl 

We note that in this study the lightest singlet scalar 
(/o) is not substantially lighter than the non-singlet (do), 
as shown in fig [3] for instance. In our previous work at 
(3 = 5.2 using the clover action with csw = 1-76 [28J, the 
singlet mass was around 50% lighter than the non-singlet 
scalar mass. 

Necco |33j has recently compared the continuum limit 
for various quantities, including the mass of the ++ glue- 
ball state, in quenched QCD using a number of improved 
gauge actions. In figure |5] we plot the masses obtained 
from the C390 and C410 data sets with the quenched 
++ glueball masses from Necco [33|. This shows that 
a lower mass is obtained with dynamical quarks for this 
improved gauge action also. 

The plot of the ++ mass as a function of the square 
of the lattice spacing (figure [3} is slightly misleading for 
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FIG. 5: FS mass from this calculation compared to the 
quenched ++ glueball mass from Necco [3^|. 



FIG. 7: Two pion threshold for the FS scalar meson as a 
function of pion mass squared. 
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FIG. 6: FS ++ mass as a function of pion mass squared. 



unqucnchcd QCD, because it does not show the quark 
mass dependence. In figureElwe plot the FS ++ mass as 
a function of the pion mass squared. We didn't attempt 
a chiral extrapolation. 

We have not included in the summary figure^ the re- 
cent results from ^| . This calculation does not see a de- 
viation of the unqucnched FS ++ masses (using glueball 
interpolating operators) from the quenched ++ glueball 
masses, but did not include simultaneous fits to the glue 
and fermionic interpolating operators. 



The dominant hadronic decay of the flavour-singlet 
q++ mesorl i s to two pions. In figure we show the 
relevance of the decay threshold. The proximity of this 
threshold may have significant impact on the spectrum, 
also it allows an investigation of the decay transition 
strength itself. 

The fo meson has an allowed S-wave decay to tttt and 
this has a low threshold so will have a significant impact 
on the analysis. In a finite spatial volume, as on a lattice, 
the two body states are discrete and one can concentrate 
on the lightest state, with the two pions at rest. A study 
of the FS sector should then involve gluonic operators 
(glueball-like) ; fermionic (qq) and two meson (tttt) to give 
the fullest coverage. We have studied the spectrum pro- 
duced by the first two types of operator above. Here we 
discuss the prospects for including the two pion channel. 

Historically the two pion channel was included in a 
quenched study of the scalar glueball by Weingarten and 
Sexton [l9| although their results were very preliminary. 
More recently two body channels have been studied |2j| 
Ej in connection with mesonic decays, such as p — > tttt 
and do — * A" A". 

In figure[7]the threshold for scalar decay to two lattice 
pions is plotted. The U355 data has lattice fo states at 
0.51(2) and 0.92(6) (in lattice units) compared to the 7T7r 
threshold on that lattice at 0.59. Thus the heavier fo is 
unstable. This implies that the impact of the two pion 
channel could potentially be strong. 

The fact that the fo and two pion states are close 
means that we may be able to compute decay widths, 
as was attempted by the GF11 group 0], using the for- 
mulation developed for the decay of the heavy 1 h hy- 
brid and p meson polli^ . 
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FIG. 8: The quark structure of the disconnected contribution 
(Q) and the triangle contribution (T) to the transition /q — » 

7T7T. 



Of course the theoretical problem of dealing with un- 
stable particles has been solved in principle by Liischcr 
over ten years ago [4^ |. In his formalism, the scattering 
phase shift can be extracted from the volume dependence 
of two-body energy levels. The mass and width of the 
resonance can then be extracted from the phase shift. 
Sufficient precision has not been available to pursue this, 
so far. Thus the more qualitative treatment of hadronic 
transitions, described above, has been used, instead. 

The lack of progress in the study of /o meson decay 
to 7T7T is due to the presence of disconnected diagrams in 
the decay- see fig. 00 Following the methods used for p 
decay, one can estimate the relative importance of this 
disconnected contribution Q compared to the connected 
triangle diagram T. The triangle contribution is rela- 
tively easy to evaluate unlike the disconnected contribu- 
tion. In fact results for the transition T of a scalar to two 
pseudoscalar mesons have been presented [2^ elsewhere. 

Using time-slice stochastic sources (as used by ref. Q) 
from 100 gauge configurations of U355, and combin- 
ing these with the volume stochastic source results used 
above, we can evaluate these contributions: see fig. [5] 
This illustrates that even though we have measured the 
disconnected contribution Q from every time-value and 
space-value on each of 100 gauge configurations, the re- 
sulting error is very large. 

Thus a thorough study of decays of FS scalar mesons 
must await a study with much larger statistics than that 
we have available. 



VI. DISCUSSION AND CONCLUSIONS 

The inclusion of the tpip ++ operators with the 
pure glue ++ operators in the variational analysis has 
produced a suppression of the masses relative to the 
quenched results at an equivalent lattice spacing (see fig- 
ure . There are two main possible interpretations of 
these results. It could be possible that the lattice errors 
(0(a 2 )) in the FS ++ masses are larger for unquenched 
QCD than for quenched QCD, because of the phase struc- 
ture of the theory. Another possibility is that in the un- 
quenched calculations the FS ++ mesons couple to the 
/ (980) or even to the / (400-1200), thus driving their 
mass way below the quenched ++ glueball mass at 1600 
MeV. We discuss these two possibilities in turn. 



2 — 
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FIG. 9: The relative contribution of the disconnected contri- 
bution (Q) (vacuum subtracted) and the triangle contribution 
(T) to the transition fa — > tvtv with all states at zero momen- 
tum for lattice U355. 



In quenched QCD the 0(a 2 ) lattice spacing errors to 
the mass of the ++ glueball are large (see the quenched 
data in figure 0J. There are lattice formalisms that add 
irrelevant operators to the lattice actions that in prin- 
ciple can remove the 0(a 2 ) errors (this is known as im- 
provement). The first systematic use of the improvement 
formalism for the ++ glucballs was by Morningstar and 
Peardon [t| . They also found large lattice spacing de- 
pendence in the mass of the ++ glueball even with their 
improved gauge actions and this phenomena is colloqui- 
ally known in lattice QCD circles as "the scalar dip". 
Further studies of the ++ masses on the lattice with 
a variety of improved gauge actions by Necco |3i| and 
Niedermayer et al. pfl |47| still found a strong lattice de- 
pendence of the results. The lattice spacing dependence 
of the ++ mass is probably the best known example of 
the break down of the simplest application of the im- 
provement and perfect action programs to reduce lattice 
spacing dependence. 

In quenched QCD the explanation of the strong lat- 
tice spacing dependence of the mass of the ++ state is 
the influence of the bulk phase transition in the adjoint 
plane 0, close to the region where calculations are 
done. The way to reduce the lattice spacing dependence 
is to add an adjoint term to the gauge action with a nega- 
tive coefficient to work in a region of the parameter space 
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away from the phase transition |3J,|50(. This matters for 
unquenched calculations because it has been conjectured 
that the clover fermion action induces an adjoint gauge 
term that will increase the lattice spacing dependence of 
the FS ++ mesons [3l|, |33, y(j . If this interpretation is 
correct, then the suppression of the FS ++ masses in 
figure 01 is purely a lattice artifact. 

There is additional evidence that suggests that the 
suppression of the masses of FS ++ mesons seen in 
this calculation is just not a lattice artifact. The ad- 
joint term is also thought to be the cause of the phase 
transition seen by JLQCD [3(| and Farchioni et al. [3?| . 
The phase structure of the unquenched lattice calcula- 
tions that use Wilson like fermion actions have been re- 
viewedby Shindler Both JLQCD [^j and Farchioni 
et al. 37] found that the phase transition was weakened 
or removed by the use of an improved gauge action, such 
as the Iwasaki action. The consistency of the suppression 
of the FS ++ masses between calculations that use the 
Wilson and Iwasaki gauge actions in figure seems to us 
good evidence that the suppression of the masses of the 
FS ++ mesons is not just a lattice artifact. 

There are other physical reasons why intuition based 
on the phase structure of quenched QCD may be not 
be a good guide for unquenched lattice calculations of 
the FS 0++ mesons. If the / (980) or / (400-1200) have 
any quark - anti-quark components then they will cou- 
ple to our interpolating operators and drive the ground 
state of the ++ channel to the 1 GeV level. Although 
there have been recent claims by Mathur et al.[Tjg that 
the ao (980) is molecular, an unquenched calculation by 
UKQCD provided evidence for the ao(980) state to be a 
quark - anti-quark state |29| . 

One of the advantages of lattice QCD calculations is 
that they can be used to understand the physical mecha- 
nisms behind the numbers. In the quenched theory the qq 
and glueball states couple to distinct ++ states. Calcu- 
lations that include very heavy dynamical quarks should 
be similar to the distinct quenched states. As the mass 
of the sea quarks is reduced the states will start to mix. 

For example: taking ao masses from ref. p9j . we get at 
U355 r m(a ) = 3.23(20). Then at this lattice spacing, 
the quenched glueball mass is around rom(GB) = 3.6. A 
mixing shift of roAE = 0.63 would move these levels to 
2.6 and 4.2 respectively, in excellent agreement with our 
two observed levels at 2.57(10) and 4.63(30) for U355. 
Furthermore, since the mixing shift is large compared to 
the initial splitting (0.37 in the above example), we would 
expect approximately maximal mixing in the observed 
spectrum - which is indeed what we find from the fit 
coefficients. Defining a mixing from 

a (0|G|l)(0|gg|2) 
a ~ <0|M|1>(0|G|2> 

where G and qq refer to the operators used to create 
the states and 1 and 2 refer to the two lightest states 
observed on the lattice. The (0|G|1) are one of the c™ 
coefficients in equation [21 For U355 with momentum 



n.n = 1, we obtain a mixing angle a = 57(10)° where 
45° would correspond to maximal mixing. 

In terms of a mass mixing matrix, this example would 
correspond to 

r m(qq) r X \ f 3.2 0.8 \ 

r X r m(GB) J \0.8 3.6 J 

so that the mixing matrix element is r$X = 0.8 and the 
resulting eigenvalues are those we observe (2.6 and 4.2) 
on the lattice. The mixing would be nearly maximal (ac- 
tually 38°) in this example. This same mixing matrix ele- 
ment (X) and glueball mass with r m(qq) = 3.56(14) [2j| 
for U350 would yield eigenstates at 2.8 and 4.4 where 
the former is in good agreement with xour lattice result, 
2.85(19), for U350. For the C410 data, at a coarser lattice 
spacing, a larger value of vqX ps 1.1 would give a split- 
ting of the r m(qq) = 3.1(1) and r m(GB) ps 3.1 values 
to agree with our lightest /o state at r mi = 2.0(1). In 
an earlier study with Nf = 2 at a coarser lattice spacing, 
an estimate of r X w 1.5 (3.65V2 0.3) was quoted J26j. 

Thus we find that in our lattice spacing and quark mass 
range, the fermionic (with Nf = 2 degenerate sea-quarks) 
and gluonic operators create states which are maximally 
mixed, with a mixing matrix element of tqX » 0.8 which 
corresponds to an energy of approximately 320 MeV. 
This value is smaller than that estimated previously us- 
ing Nf = 2 sea-quarks with a coarser lattice spacing and 
than that we estimate from C410 with a similar coarse 
lattice spacing. 

If the mixing matrix element X were to remain con- 
stant as the quark mass is reduced, then the estimate |29l ] 
of the ao mass (for Nf — 2) is 1 GeV while the glueball 
is around 1.6 GeV. These input masses would then be 
mixed to 0.86 and 1.74 GeV. This gives some indica- 
tion of the magnitude of mixing effects which could be 
present in the experimental spectrum and require calcu- 
lation with similar small lattice spacings. 

In the experimental spectrum, however, there are ex- 
pected to be significant effects arising from the ss scalar 
meson, which we have neglected here. This ss state 
should mix with gluonic operators and one will need a 
3x3 model to incorporate this adequately. Lattice QCD 
is able to include 2+1 flavours of sea-quark and this will 
be appropriate for a fuller confrontation with experiment. 

In quenched QCD the first calculation Q that com- 
puted the mass of the ++ glueball in the continuum 
limit used a lattice spacings down to 0.05 fm to have a 
reliable continuum extrapolation. Unquenched calcula- 
tions that use the Wilson gauge action may also have 
similar scaling violations. 

As we have noted, there appears to be a significant 
lattice spacing dependence in the mixing. To quantify 
the amount of glueball to qq mixing in /o mesons will 
require a dynamical fermion calculation at even finer lat- 
tice spacings such as 0.05 fm. For Wilson type fermions, 
recent algorithmic advances [5^, IBH. l54| means that this 
may now be just attainable on the current generation of 
machines. 
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The eventual definitive study of flavour singlet scalar 
mesons will need fine lattice spacing, dynamical simu- 
lations with light quarks (plus a strange quark in the 
sea) and large statistics to enable the disconnected dia- 
grams to be evaluated accurately. The methods we have 
used with Nf = 2 degenerate sea-quarks of mass down 
to about 50% of the strange quark mass show the way 
forward. We find substantial mixing of glueballs with qq 
states in the fo spectrum and we expect that feature to 
survive in a future study. 

Although, given the qualifications we mention above 
about this calculation, the message we get from the sum- 
mary of glueball masses in figure 0] is that the effect of 
unquenched fermions is to drive the mass of the FS state 
from the 1600 MeV of the quenched glueball towards 1 
GeV. This calculation does not favour the weak mixing 
of glue and ipip operators proposed by Weingarten and 
Lee and other groups [5^, . Looking at our data it is 



not clear that a mixing scheme based on only the states 
/o(1370), /o(1500), and / (1710) is complete enough to 
determine the fate of the quenched glueball. 
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